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Abstract 

We study linear stability of planar travelling waves for a scalar reaction- 
diffusion equation with non-linear anisotropic diffusion. The mathemati- 
cal model is derived from the full thermo-hydrodynamical model describ- 
ing the process of Inertial Confinement Fusion. We show that solutions 
of the Cauchy problem with physically relevant initial data become pla- 
nar exponentially fast with rate s(e' , k) > 0, where e' = ^ min <C 1 is a 
small temperature ratio and > 1 the transversal wrinkling wavenum- 
ber of perturbations. We rigorously recover in some particular limit 
(s',k) — > (0, +oo) a dispersion relation s(e',k) ~ "fok a previously com- 
puted heuristically and numerically in some physical models of Inertial 
Confinement Fusion. 

1 Introduction 

The aim of Inertial Confinement Fusion - ICF in short - is to concentrate 
extreme temperature and pressure conditions in plasma to a single point in 
order to trigger a nuclear fusion reaction. The corresponding models couple 
thermal reaction-diffusion equations with hydrodynamical effects. For such 
extreme temperature regimes (T ~ lO 7 ^) the heat conductivity coefficient A 
in the usual diffusion term V ■ (AVT) cannot be considered as constant, but 
must be taken of the form 

A = A(T) = T m_1 

for some fixed conductivity exponent m > 1 (the value m = 7/2 is often 
considered in ICF, corresponding to the Spitzer electronic heat conductiv- 
ity). Scaling at the hot side we always take in the sequel temperature values 
e' = T min < T < T max = 1, and consider the limit of small temperature ratio 
e' < 1 (which is indeed true in ICF, typically iH 3 - ~ 10~ 3 ). An important 

v J max ' 

consequence is here the strong variation of the diffusion length scale: since 
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temperature varies of several orders of magnitude so does the thermal con- 
ductivity A = T m_1 , and the effective diffusion length scales are thus very 
different in the cold and hot regions. The physically relevant wavenumbers 
are therefore 1 <C k <C l/(£') m_1 , corresponding to wavelengths h between 
the shortest and longest diffusive length scales. The second important feature 
of 1CF is the propagation of spherical temperature waves, which concentrate 
some laser driven input of energy to the center of a fuel target. These waves 
have a very particular spatial structure consisting in three separate zones: in 
the most inner one (the target) an initial amount of unburned fuel is at rest 
T fn e', in the outer one the fuel has been transformed into hot plasma by the 
laser T« 1, and in the intermediate region temperature evolves between e' 
and 1 according to a specific power law. 

When the temperature ratio e' is small a very thin boundary layer, called 
ablation front, separates in the wave the cold region from the intermediate one. 
In the theoretical situation the geometry is purely radial and the ablation front 
should be a sphere shrinking in time towards its center. This ideal scenario is 
however unstable due to Ray leigh- Taylor and Darrieus-Landau effects, which 
are purely hydrodynamical and concentrate at the ablation front: should the 
latter slightly deviate from its ideal spherical geometry, wrinkling may amplify, 
centripetal waves may stop propagating, and the desired fusion reaction may 
not be triggered. A so-called transversal mass ablation effect fortunately tends 
to stabilize the system. This phenomenon can be roughly explained as follows: 
whereas the diffusion in the radial direction is responsible for the particular 
structure of the wave, the orthoradial diffusion across the wrinkled front simply 
tends to reduce its unstable transversal oscillations. 

A self-consistent linear analysis of the full thermo-hydrodynamical model 
[llj, performed in the relevant frequency regime for large conductivity expo- 
nents m +oo, showed that the transversal mass ablation should lead to a 
stabilizing contribution 

s(e', k) ~ z/fc 1- ™ 31 , v = est > 

in the full dispersion relation (s > corresponding to an exponential damp- 
ing rate). It was suggested in [8] that this stabilization mechanism can ac- 
tually be investigated looking at a much simpler model, namely the purely 
diffusive relaxation of thermal waves, free from hydrodynamical effects. A 
first rigorous study of an approximated non-Galilean model [6j for finite con- 
ductivity exponents led to a very similar dispersion relation in the regime 
1 <C k 2 <C l/{£ , ) m ~ 1 , except for a logarithmic correction. In this paper we de- 
rive a purely thermal model suitably approximating the full one, and rigorously 
recover the anticipated self-consistent dispersion relation s ~ z^A; 1- ™^ 1 for fi- 
nite conductivity exponents m > 3 encompassing the physical case m = 7/2 
(in some frequency regime discussed later and relating k 3> 1 to s' <C 1). 
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2 Model and contents 



The spherical geometry being difficult to investigate we perform our study for 
the following planar configuration: x G R will denote the infinite longitudi- 
nal direction of propagation (mimicking the radial depth of penetration in the 
target for the real setting) and y £ R the transversal one. This planar approx- 
imation is legitimate in ICF since the relevant wavelengths are much smaller 
than the radius of the target. Neglecting the transversal velocity V = (V, 0) 
we consider now the following 2 dimensional thermo-hydrodynamical model 
for (t,x,y) e M 3 

pd t T-W ■(XVT) + pVd x T = f(T), (2.1a) 

a tP + d x (j>v) = o, (2.1b) 

pT = 1, (2.1c) 
T(t,-oo,y) = e' T(t, +oo, y) = 1. (2.1d) 

Here T is the reduced temperature of the plasma scaled at the hot side 
(Tmin = e' « 1 = T max ), p its density, V G R its longitudinal velocity, 
A = A(T) = T m ~ l the reduced non-linear heat conductivity, and f(T) a non- 
linear reaction term of ignition type modelling input of energy by the laser 



in the outer region. Equations (2.1a)-(2.1b) correspond to conservation of en- 



ergy and mass, and (2.1c) is a quasi-isobaric approximation p w est which is 
common in ICF. Since we neglected the transversal velocity conservation of 
momentum is not required here to close the system, which is however Galilean 



invariant. Boundary conditions (2. Id) simply state that the medium is at rest 
at x = — oo and totally combusted at x = +oo. Temperature waves propagate 
now in the negative x direction (the hot medium invades the cold one), and the 
transversal wrinkling of the front will be later taken into account considering 
a periodicity condition in the y direction. 



Considering ( |2.1b ) as a Schwarz condition for crossed derivatives we may 



define X(t,x,y) such that d t X = —pV and d x X = p. More explicitly, this 
new Lagrangian coordinate X is mass-weighted with respect to the Eulerian 

one x as ^ 

X{t,x,y) := / p{t,z,y)dz- I pV(r,0,y)dr. (2.2) 
Jo Jo 



One easily checks that for physical p > 1/T > 1 (2.2) defines a proper change 



of coordinates (t, x, y) — > (t, X, y) such that X(t, ±oo, y) = ±oo, and inverting 
the Jacobian yields the usual characteristics || = V(t, x, y) for Lagrangian 
particles. The ideal spherical geometry corresponds here to planar longitudinal 
waves d y = 0: from now on we use the 

Lagrangian approximation: d y X = 0. (2-3) 

Since coupling of the conservation of mass and energy occurs precisely through 
this d y X term this is of course a drastic approximation, but however more 
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subtle than neglecting each and every transversal variation d y = directly in 
the equations: transversal diffusion will still be effective and indeed induce the 
stabilization mechanism as explained in the introduction. Throughout the rest 
of the paper we will work in these Lagrangian coordinates and write x instead 
of X with a clear abuse of notations, and in the Lagrangian approximation 



(2.3) conservation of energy (2.1a) is uncoupled from the hydrodynamics as 



pd t T - pd x (Xpd x T) - d y (Xd y T) = f(T). 



(2.4) 



Since p = 1/T and A = T rn (2.4) is a purely scalar react ion- diffusion equa- 



tion, in which it is worth noting the clear difference between longitudinal and 
transversal diffusions. 

As previously mentioned, temperature in the wave profile evolves according 
to a specific power law. More precisely, v = T m ~ l behaves linearly with respect 
to the Eulerian longitudinal coordinate in this intermediate region: this is very 
similar to the Porous Media Equation where the pressure unknown evolves 
linearly at least in the vicinity of the free boundary, see e.g. [5]. However, the 
Lagrangian coordinates are here mass-weighted and stretched with respect to 
the Eulerian ones according to (2.2): the suitable variable turns to be 



T 



m-2 



which we will show to behave linearly with respect to the Lagrangian coordi- 
nate at least in the intermediate region (existence and qualitative properties 
of the wave solution will be investigated in Section |3|. Since we are concerned 
with wave propagation it will also be more convenient to work in the moving 



frame x + ct, in which conservation of energy (2.4) and boundary conditions 
(2. Id) read in terms of the unknown p = T m ~~ 



d t p + cd x p - ( pd xx p + 



1 



p(t, -oo, y) = e, 



m — 2 
p(t,+oo,y) 



{d x pf 



as 

m 

P™-*dyyP + 

l\m-2 



GOO 



£ : 



(ej 



< 1. 



(2.5a) 
(2.5b) 



The right-hand side G(p) := (to — 2)/i/( y u™ 3 2) corresponds to the initial 
reaction term f(T) modelling the laser, and is again of ignition type: we 
will always assume in the following that G is smooth and that there exists a 
threshold 9 G]0, 1[ such that 



tie [0,6] 
pe]9,l\ 
p = i 



G = 0, 
G > 0, 

G = and G'(l) < 0. 



In order to investigate linear stability with respect to wrinkled pertur- 
bations U(t,x,y) we linearize in the moving framce at the planar travelling 
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wave /j,o(t,x,y) = p{x), with an additional 2 r njk periodicity condition in the y 
direction accounting for the ablation front transversal wrinkling: 

d t U - (pd xx U + p^d yy u) + (c - d x U - p"U = ^{p)U, (2.6a) 

U{t, ±00, y) = 0, (2.6b) 
U(t,x,y + 2ir/k) = U(t,x,y). (2.6c) 

Note again the difference between the longitudinal diffusion pd xx and the 

m 1 1 

transversal one pm-2 Q yy i n (|2.6a|). Expanding now in transversal Fourier modes 



U(t,x,y) = u(x)e st+lk y leads to 

— pu" + I c — ) u + \k 2 p"^ — p" — (p) ] u = su, (2.7a) 

\ m — 2 J \ dfj, J 

u(±oo) = 0, (2.7b) 

where ' = d/dx. 

Let us stress at this point that the wave speed c = c £ > and profile 
p = p £ (x) > are uniquely determined by the small parameter e > (this 
should be no surprise since the reaction term is of ignition type, see Section [3]). 
We will see that the planar wave satisfies lim p(x) = e with exponential rate 

X— > — 00 

r = c £ /e: when £->0we have r — » +00 and the wave degenerates into a free 
boundary solution, p £ — > p and po(x) = for x sufficiently negative. This is 
illustrated in Figure [T] and is again very similar to the Porous Media Equation 
scenario. Physically relevant perturbations should therefore not disturb the 
reference wave ahead of the front, and as a consequence we only investigate 
perturbations with maximal decay u(—oo) = 0. 

For fixed temperature ratio e > and wavenumber k > only some 
particular values of s = s(e,k) allow the connection u(— 00) = u(+oo) = 
with maximal decay on the left side, and since we set s to be the damping 
rate we also want to compute the smallest such s (corresponding to optimal 



stability). Since the wave profile will satisfy p(x) > e > 0, problem (2.7) is 
uniformly elliptic on the line x G M: finding the minimal s is clearly a principal 
eigenvalue problem, and signed solutions will play an important role as usually 
in classical elliptic theory. 

Let us recall that the relevant frequency regime in ICF is 1 k 
l/i^e')" 1 ^ 1 = l/em-2; in the whole paper we will always assume the condition 

1<A;<^^, (2.8) 

£m-2 

which is slightly more restrictive for wavenumbers if m > 2. For some large 
k > and small 77 > let us define the frequency (open) set 

U:=\e,k>0: k > k and k < — — \. (2.9) 
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Taking (e,k) — > (0,+oo) in U clearly implies regime (2.8). The parameter 



7] > is needed for technical reasons, but may be chosen as small as desired 



so that there is no real loss of information when restricting (2.8) to (e, k) G U. 
Our main result is the following: 

Theorem 1. Fix a reaction term G of ignition type and a conductivity expo- 
nent m > 3: for any n > there exists fco > such that, for any (e, k) G U 



defined in (2.9), there exists a principal eigenvalue s = s(e,k) > such that 
problem (2.7) has a positive solution u(x) with maximal decay when x — > — oo. 
Moreover s(e, k) > is the smallest eigenvalue, there exists no other principal 
eigenvalue, and the associated eigenspace is one- dimensional. Finally there 
exists 70 = 7o(^i, G) > such that 

s{e,k) ~7o/c 1- ^ (2.10) 

when (e, k) — > (0, +oo) in U. 

This means that the wave is linearly stable with respect to wrinkled pertur- 
bations (with maximal decay) and that any solution of the non-linear Cauchy 
problem with initial datum Uq(x) decaying fast enough to the left will become 
planar exponentially fast with rate at least s(e, k) > 0. Note that our result 
does not tell anything regarding perturbations in the longitudinal direction: 
invariance of the wave solutions under x shift results as usual in a non-trivial 
element dp/dx in the kernel of the linearized operator, i-e in a zero eigenvalue. 
For fixed e > non-linear techniques can be used to establish stability of the 
wave with respect to almost any such perturbation (see e.g.jHl HO]), but one 
has then to face the degeneracy into free boundary when e — > 0. The main 



interest of our result is that it is uniform in e — > 0, see (2.10). 

The paper is organized as follows: in Section [3] we establish existence of 
the planar travelling wave and investigate its linear behaviour as well as the 
degeneracy into free boundary of the ablation front when e — > + . Section [4] 
is devoted to the construction of solutions with maximal decay at — oo in the 
cold region p = O(e) for s = ^k 1 ' m - 1 and fixed 7. This is done expanding the 
solution u = uq + eu\ + eu2 in the suitable scale £ = | , and we also compute 
(ti, ^) at the exit of the boundary layer. Section [5] is a formal limit (e, k) — > 

(0, +00) performed in the linear region ^ fa est > at scale ( = k 1 ^ m -^x. Due 
to the presence of the free boundary when e — > the resulting limiting problem 
degenerates at ( = 0, but will nonetheless yield the asymptotic coefficient 70 



in (2.10) as the principal eigenvalue of some singular Sturm-Liouville problem 
on the half line ( > 0. In Section [6] we rigorously justify this limit by proving 
that, when (e, k) — > (0, +00) in U, the real physical setting e > 0, k < +00 
automatically matches the formal asymptotic problem e = 0, k = +00 (the 
matching occurring in the ablation front). This last section also contains the 
proof of Theorem [TJ 
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3 Planar travelling wave and boundary layer 

In this section we study the aforementioned travelling wave at which the prob- 
lem is linearized. Plugging /i(t,x,y) = p(x) into (2.5a) leads to 

-pp"+(c-^-)p' = G(p), p(-oo)=e, p(+oo) = l, (3.1) 
\ to — 2 j 

and in order to eliminate invariance under x shift we impose the additional 
pinning condition 

p(0) = 9. (3.2) 

Proposition 3.1. Let e e]O,0[: 



1. There exists a unique speed c = c £ > such that problem (3.1)-(3.2) 
admits a solution p = p £ (x). This solution is unique and satisfies < 
p' e < (to - 2)c E . 

2. When e — > + we have that c £ — > Cq > and p £ {.) — > Po(-) uniformly on 
K ; where po solves 

xG ]-° '-{^)^] : Po(^) = 

: PoW = ^ + (m - 2)c x 

-PoP , o , + (co-^)p^ = G(p ) 
Po(0) = ^, Po(+oo) = 1 



X 



x elo, +oo] 



This is a particular case of more general well-known results in reaction- 
diffusion theory, see e.g. |21 H]: the proof relies on simple ODE techniques and 
will only be sketched here for the sake of completeness. The degeneracy into 
free-boundary p £ (x) — » po(x) is portrayed in Figure [TJ 

Proof. Taking advantage of the invariance of ( |3.1[ ) under translations we may 
use the Sliding Method from [3] and see that any solution must be increasing 
in x, and therefore that e < p < 1. Setting p' = U{p) > 0, (3.1) leads to 



pe]e,l[: 



U(p) 
m — 2 



C/(p) = G(p) 



with C/(e) = p'(— oo) = and U(l) = p'(+oo) = 0. For p < 9 the reaction 
term can be omitted and this can be explicitly integrated as 



pe]e,9]: U(p) = (m-2)c 
and the pinning condition p(0) = 9 yields 



1 



(3.3) 



p'(0) = U(9) = (m - 2)c 



1 - 
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For fixed c > and given a > the Cauchy problem 



-pp 



x>0 



p(0) 

p'(o) 



G(p) 





a: 



is independent of e, and shooting with respect to a there exists a unique 
a = «o(c) > such that the corresponding solution satisfies p(+oo) = 1. 
The solution on x < therefore matches the one on x > if and only if 

(m — 2)c ^1 — (|) m - 2 j = p'(Q) = «o(c)- One then shows that this fixed point 

equation in c has a unique solution c = c e («o is decreasing in c) associated 
with a unique profile p = p £ , and convergence of c e — >• Co and p e — >• po when 
£ \ follows by monotonicity. □ 



x 



We establish now precise asymptotes on p £ — > p when £ — > 0. 



e£ and 



Scaling 



(3.3) reads 



dq 



dp 
dx 



q 



m-2)c[ l-—^- 

qm~ 

OO) = 1 



This shows that p £ presents a boundary layer of thickness e in which p — 0(e) 
(this is the ablation front), and p' varies from p'(— oo) = when p = e (thus 
g = 1) to (m — 2)c[l + o(l)] > when p ^> £ (thus g ^> 1) as pictured in 
Figure [1} 





1 


Po(x) 













. /• 


e 














O(e) 



Figure 1: monotonic convergence p e — > po and £-boundary layer. 



From now on we choose the origin at the slope discontinuity by setting 

q 

x e-=7 ^v— > Po(aJ«)=0, (3-4) 

(m - 2)c 

and we are out of the boundary layer as soon as x > £. To the right p e is 
linear when e — > 0, p'Jx) ~ (m — 2)co and p e (x) ~ (m — 2)cox. We will 

e— >0 e— >0 
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refer in the following to the set where e <C p < 6 indistinctly as the "hot zone" 
or "linear zone", and call "cold zone" the set where p = 0(e) (Figure [l] should 
make this terminology self-explanatory). In addition, plugging ( |3.3[ ) into ( |3.1 ) 
yields an expression of 

p ,^LLuM\ u() ^ CE ^uM { ey (3 . 5) 

p £ V m - 2 J p £ ^ 



for any e > 0. As a consequence of (3.3)-(3.5) we obtain the following asymp- 
totes in the linear zone when e — > + : 

p e (x) ~ (m — 2)cox, (3.6a) 
p' £ (x)~(m-2)co, (3.6b) 

with A := c ° j > 0. These only hold far enough from the boundary 

[(m-2)co] m " 2 

layer a; ^> e and as long as p £ < 9 (otherwise the reaction term must be taken 
into accounted). 



Finally, asymptotic analysis of (3.1) at x — — oo (where p = e and p' = 0) 



easily yields for fixed e > the exponential decay 

p' £ (x) ~ e c ^/ £ , (3.7) 



and \p £ (x) — e\ = O (e°s x ). This will be useful in the next section when 

— oo 

building the maximal decay solution in the cold zone. Note again how the 
rate blowup c £ /e ~ c /e — > +oo translates the degeneracy into free-boundary. 



4 Cold zone and asymptotic expansion 



In this section we construct the solution « of (2.7a) with maximal decay in 
the cold zone in the form of an asymptotic expansion u = Uq + eu\ + eu2, with 
u 2 U\ in some sense. For the sake of simplicity we omit the subscripts and 
write here c = c e and p = p e . Let us anticipate that perturbations will have a 
spatial structure resembling the one of the wave, including a boundary layer 
of thickness 0(e): we therefore scale again as 

x = q(0 = -p(et), «(0 = u(e£)- 

The slope of the reference wave solution jumps inside the boundary layer from 
dp/dx = dq/d^ ~ to dp/dx = dq/d^ (m — 2)c > 0, and this transition 
is steeper and steeper when e — > + : we expect a singularity of the second 
derivative somewhere, which is of course consistent with the slope discontinuity 
in the asymptotic profile po(x). In order to later neglect p"(x) in the linear 
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zone we will have to push the exit of the boundary layer far enough so that our 



asymptotic expansion encompasses this singularity. Asymptote (3.6c) shows 

i 

that in order to do so we need setting 



x £ :=e l - a } &:=x e e, a := 1 + M 4.1 

777 — 1 V 2/ 

1 

for some 77 > 0, it is easy to see that e m ~ 2 < i e < 1 holds if rj > 0. As a 
consequence we may safely neglect p"(x) for x ^> x £ . Here rj > is exactly 
the one in Theorem [T] and arbitrarily small, but its occurrence in the above 
definition of a is purely technical. In agreement with Theorem [T] we anticipate 
that the relevant values of s should be of order k 1 ^™^ 3> 1: through this 
whole section one should think of s as of for some fixed constant 7 

of order unity and independent of e, k, which will be adjusted later. We also 



recall that we are interested in (e, k) — > (0, +00) in the double limit (2.9). 
Since we set x e 3> e out of the boundary layer we have p(x e ) ~ (m — 

2)cqX e 1, and we may therefore omit the reaction term (of ignition type) on 



the cold interval I =] — 00, x e ]. Recasting (2.7a) in £ coordinates as 



Lv = ehv 

.f (,. 2'/' 



L--=-Q^+{c-^r 2 )±- q », (4 . 2) 
h = I s — k z e m - 2 q m - 2 ) , 



we shall seek below solutions in the form v = Vq + ev\ + evi with v 2 v\ 
in some sense. We will refer to L as the principal operator, which is also the 
linearized operator for non- wrinkled perturbations k — 0. 

4.1 Maximal decay and principal operator 

Let us recall that we investigate maximal decay perturbations u(— 00) = 0: 



the asymptotic equation associated with (2.7a) for x = —00 yields two char- 
acteristic exponential rates 



c ± \ c 2 + 4e(k 2 E™- 2 — s) 



For s of order k m - 1 regime (2.9) implies that < s(s — k 2 e m - 2 ) 1, and 



(4.3) consequently yields 



c c 
< r~ < r + ~ -, r + < -. (4.4) 

e e 



Maximal decay obviously corresponds here to u (x) = o(e r+x ^j. 
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In order to work in some fixed functional setting (independent of e) let us 
define the following weighted spaces on / =] — oo, £ e ] 



r , <n K = {/eft, /»eC 6 (/)}, 

= i i k and B - = {/ G c "> v ^ fc : / w G B ° ) * 

1 4 <o = {/ G /(&) = 0} 

(recall that c £ — > Cq > 0). Here C\ denotes the space of functions with contin- 
uous bounded first k derivatives, and the B w J s are equipped with their usual 



Banach norm. On one hand (4.4) shows that any maximal decay solution 
must behave as e er+ ^ when £ — > — oo with er + ~ cq > Cq/2, and therefore 

e,k 

belongs to B^. On the other hand any non maximal decay solution behaves 
as t>(£) ~ e £r * with < er~ <C Cq/2: such solutions cannot belong to B^, 

-oo e,k 

and as a consequence it is legitimate to look for maximal decay solutions only 
inR ,. 



w ' 



Lemma 4.1. The principal operator L : B^ — > B^ is continuously invertible 
and 

\\L- 1 \\<Ce- a . (4.5) 



Proof. Defining $ =: J - it is easy to check that Duhamel's formula 



/(0 := 




exp[$(r/) - $(z)] 



drj dz 



(4.6) 



yields the unique solution in B^ of Lf = g for any given g G B w . The estimate 



for L 1 then follows from (3.6)-(3.7) expressed in £ = x/e coordinates and the 
study of $ for £ — oo and 1 <C £ < £ e = e _a combined with (4.6). □ 



4.2 Asymptotic expansion and frequency regime 

Vq + ev\ + et> 2 with f 2 <C i>j in some sense leads to solving 



Expanding v 

Lv = ehv as 



Lvq = 0, 
Lv\ = hvo, 
[L — eh]v 2 = ehv\. 



(4.7a) 
(4.7b) 
(4.7c) 



In order normalize perturbations v(£ e ) = 1 we will require in addition that 



Resolution of (4.7a): the reference wave solution is as usual determined 



up to shifts and L[q'] = 0, thus yielding a suitable candidate v = q' 
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for the leading order. This is confirmed by (3.7), showing that q'(£) 



p'(s£) ~ e <C e er 5 decays fast enough. Normalizing finally leads to 

— oo — oo 



■ 1 / m m \ 

Resolution of (4.7b): remark that h = Is — k e m ~ 2 q m - 2 ) G L°°(I) im- 



plies /it>o G Lemma 4.1 then properly defines 

Vi := L _1 [/iUo] G . 



Resolution of (4.7c): taking advantage of the linear behaviour of g(£) for 
large £ and regime (2.9) with s = 'yk 1 ^™^ 1 (for fixed 7 of order unity), 



we see that h = s — k e m - 2 q m ~ 2 ~ s uniformly on / =] — 00, £ e j: we may 
therefore consider the restriction of L _1 to. The (continuous) restriction 
of L^ 1 to the subspace -B^ therefore satisfies 



\eL~ l [h- 



l£«o) 



<e\\L- l \\. 



< Ce^s < 1 (4. 

£,k 



according to Lemma 4.ll and M := Id — eL 1 h G C(B%, n ) is hence close 



to Identity thus invertible. This finally allows to solve ( j4 . 7c ) as 
v 2 :=M- 1 (L- 1 (e^ 1 ))G J Bi . 



Regime (2.9) also ensures that 



= Vq + ev\ + evi is really an asymptotic 
expansion in B^, in the sense that ||eu 2 ||sa "C ||£t>i|| B 2 <C ||^o||bo = C(l). 
More precisely, one easily shows using (4.8) and M m Id that 



vi = L 1 (hv ) 
v 2 ~ L^iehvi) 



\evi\\ B o < Ce 1 a s| |^o| I ro < ||vo||fl2 



k^H^o < Ce 



1—a 



(4.9) 



4.3 Exit boundary conditions 

We compute in this section (v,v') at the exit £ = £ e of the cold zone. Since we 
normalized u(£ £ ) = 1 it is enough to compute v', and we estimate separately 

• Since we set the exit far enough out of the boundary layer, asymptotes 



(3.6) hold at x £ = e 



1 a 3> e and therefore e€^(x E ) ~ 



Ae 



o 1+: 



terms of £ = x/e this corresponds to 



In 



(4.10) 
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i L°° (I) 

Let us recall that Lv\ = Hvq = h -^k ; and h ~ s, so that t>i is close 
in B^q to the solution of Lv = s q , 9 ^ . Using (4.6) it is easy to compute 
explicitly L- 1 ^) = (& - hence 



(4.11; 



Since the -B^, topology controls the first two derivatives we obtain 



d 



e,k c £ q'(Q d£ I 



£=6 



s 

c F 



(4.12) 



Recalling that M — (id - eL^h) w Id in C(Bl fl ), ( ]47c) and jOI| ) 
show that 



W2 = ( id-eL- 1 ^ )- 1 ^- 1 ^ « ^T^y^" 1 [(6. - 0/ 



Maid 



in -B^q. Duhamel Formula (4.6) gives an explicit integral formula for 
[(£ E — £)q']: differentiating with respect to £, evaluating at £ = £ £ 
and taking advantage of the asymptotic behaviour of $ (for £ — > — oo 
and 1 <C £ < £ e ), one finally estimates 



■/ 2 (&) = O (e 1 -* 2 ) . 



(4.13) 



We claim now that the exit slope v'(£ e ) is dominated by ev[(£ s ) in the double 
limit (e, k) -» (0, +oo): 

• From (4.10) and ( |4.12 ) we see that |i>d(£ e )| <C |efi(£ £ )| holds as soon as 
5 a ( 1+ m-i) <g; £S) which is true with s = 7/c 1- ™- 1 ; e <^ \ <^ k and by 
definition of a = a=f (l + f ). 

m— 1 V 2/ 

• Similarly from (4.12) and ( |4.13 ) we see that |ei>2(£e)| "C |£^i(£e)| holds 
in the regime (2.9) with our choice of a. 



Let us recall that s is here of order k 1 m ~ 1 ^> 1 as anticipated from Theorem[l 
the new parameter 

s 

a := — 

(m — 2)c £ k ~m-i 

should therefore take values of order unity, and the exit conditions are finally 
summarized in terms of this new parameter a by 



— {m — 2)aek m - 1 



(4.14) 



For technical reasons we will also need 
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Proposition 4.1. For fixed (e,k) the quantity v' '(£ e ) is continuously differen- 
tiable with respect to a and 



d_ 



[v'(Q\ = -(m- 2)6^-—^ + r(e, k, a) 



(4.15) 



holds with r(s,k,a) = o^ek 1 m-ij when (e,k) -^4 (0, +oo) and locally uni- 
formly in a. 



Remark that (4.15) is indeed consistent with a formal differentiation of 



(4.14) with respect to a. 

Proof. Regularity with respect to a is a classical consequence of the linear de- 
pendence on a of (4.2). Writing z = estimate (4.15) is just a computation 

for *j|(£ e ) = —(m — 2)ek 1 ~ m - 1 + ... with z = Zq + ez\ + ez<i- Exactly as for 

the slope ^ the order one determines here the dependence on a at the exit, 

namely %(£e) ~ £-Jr(£ e ). The technical computations are omitted here in 



e,k 



order to keep this paper in a reasonable length. 



□ 



5 Hot zone and asymptotic problem 

In this section we take a formal limit e = 0, k = +oo in some scaled ( coordi- 
nates corresponding to a suitable (infinite) zoom in the linear zone. We obtain 
an asymptotic eigenvalue problem on the half line ( > (( = corresponds 
to the exit of the boundary layer after zooming out), and this will relate the 
asymptotic coefficient 7 in Theorem [T] to some principal eigenvalue o"o of the 
limiting problem. The main issue in this section is precisely to compute this <7o. 

We anticipated that s should be of order A; 1- ™^ 1 , which also turns to be 



the suitable length-scale to investigate the linear region: scaling (2.7a) as 
( = k l ~ a 



(m-2)ck 



i 

m — 1 



k m ~ 1 



v{() = u ( _ c t ) q{() = k 1 m-i p ( 

\ k m-l / \ 

for k < oo leads to 

— qv + I c I v — q v = I [m — 2)ca — q™- 2 + 



(5.1) 



1 . 
m — 1 



m — ZJ V k 

(5.2) 

with q = q E> k(C) an( i where p = p £ and c = c £ depend only on e. This scaling 
is again Lipschitz g = and the parameter cr takes values of order unity. 

We recall from Section [3] that in the limit e — » we have c e — > Cq > 0, and 
the asymptotic wave profile p = lim p £ is exactly linear for x G [0, ( m _ 2 ) Co ] ■ 
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In C coordinates the exit of the linear zone x = Xg corresponds to 



Ok 1 '—! 

Ck ■= 7 ^ -> +00, 

[m — 2)Co 

while the exit of the boundary layer x = x £ is now 

( £ = k 1 ~^x £ = k^^e 1 ^ — > 

e,k 



in the regime (2.9). The wave profiles behave as q(() w yo(C) = (m — 2)cq( 

dp 

11 1"/ 1 /■ ii" ^ 1 /ii" 1 1 \ i" I /v-i f tit/1 m m r v i / w v 1 / w 1 1" 

da; 2 



in the linear zone £ G [£ e , Cfc], which grows to ]0, +00 [. Moreover since ^(po) 
and are of order 1 for x G]0, +00 [, we may neglect 



jE^ m— 1 ^ m — l m—1 k 

at least formally when e = and = +00 on ( G]0,+oo[. As a matter of 
fact a singularity of appears inside the boundary layer due to the slope 
discontinuity of po, but this will be carefully examined in Section [6] (this is 
precisely the reason why we pushed the exit of the cold x = x £ zone far 
enough) . 

Formally taking e = 0, k = +00, q(Q = 90(C) = {m — 2)c C, Ce — and 



Cfc = +00 in (5.2) we obtain the following asymptotic problem: 

v' 

CG]0,+oo[, -Cv" - + bv = av, (5.3) 



m 

with 



2 



b(£) = BC^, B := [(m - 2)c ]^ > 0. 
Since we were originally interested in perturbations u(x) vanishing at infinity 



we also require in (5.3) that t>(+oo) = 0, which is again a formal but suitable 
condition. 



Remark 5.1. For fixed a and( = +00 equation (5.3) reads —v"+B( m - 2 v = 0. 
One easily sees that there exists a one dimensional space of stable solutions 
v(+oo) = 0, while any other solution must blow v(+oo) = +00 (both at least 
exponentially). 

Recalling now that we also look for signed perturbations (as principal eigen- 
functions), the parameter o clearly appears here as a principal eigenvalue to 



be determined for the problem (5.3) with associated zero boundary condition 
at infinity. Two main difficulties arise here compared to usual Sturm-Liouville 
theory: unboundedness C G IR + and singularity at C = 0. As a matter of fact 
we shall seek solutions which are C 1 at C = 0, and one should think of this 
regularity as a boundary condition. 
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Throughout the rest of Section [5] we will consider e = 0, k = +00, and the 
only relevant dependence will be the one on a. In order to keep our notations 
as light as possible we will write p = po(x) = (m—2)cox, q = <Zo(C) = ( m ~2)coC 
and c = Cq. The main result in this section is: 

Theorem 2. There exists a principal eigenvalue o"o > such that problem 



(5.3) has a positive solution vq G C 1 ([0, +oo[)flC 2 (]0, +00 [) satisfying boundary 
conditions Vq(0) = 1 and t>(+oo) = 0. This principal eigenfunction satisfies in 
addition v' (0) = — (m — 2)o"o and v' < 0. 

In this statement the C 1 regularity at £ = and v(+oo) = should be 
seen as a boundary condition, in contrast with Vq(0) = 1 which is simply a 



normalization. Note that according to scaling (5.1) o"o of order unity means s 
of order fc 1- ™^ 1 : this will yield the asymptotic coefficient 70 = (m — 2)cqOq in 
Theorem [T] 

5.1 The asymptotic principal eigenvalue 

In this section we first establish some technical results and then prove Theorem 
[2j We start by studying the singularity at ( = 0. 

Proposition 5.1. For fixed d6R there exists a unique solution v a £ C 2 (]0, +oo[)fl 



C 1 ([0, +00 [) of (5.3) such that v a (0) = 1. In addition this solution satisfies 
v' a (0) = — (m — 2) a, and the mapping a (->• v a (.) is C 1 from R into C 1 ([0, Co]) 
for any fixed > . 

Proof. Remarking that the exponent in the zero-th order coefficient b(() = 

Til 1 ■ 

i s > i ; solutions of (5.3) "should not see" this coefficient in the 



neighbourhood of £ = and therefore behave as 



w 

(w" + + aw = 0. (5.4) 

m — 2 



It is easy to obtain a regular solution w G C 2 (]0, +00 [) D C 1 ([0, +oo[) of (5.4) 
satisfying w(0) = 1 and w'(0) = — (m — 2)o~ in the form of a power series. 
This solution is unique (see [7] Theorem 6.1 p. 169, first kind singularities) and 
classically regular with respect to a. 



We may now look for solutions of (5.3) in the form v — w + h: h should 



obviously solve a non-homogeneous equation involving w, and the limit £ — > 
also yields h(0) = h'(0) = (remember that w(0) = 1 and w'(0) = —(m—2)a). 
The proof consists in two steps: a classical Banach fixed point on the integral 
formulation for h first yields a solution on [0, Q for fixed small ( > 0, thus 
stepping away from the singularity. This solution then extends to the right 
[(, +00 [, where the ODE is now regular. Uniqueness is obtained as for w since 
the singularity at ( = is of the first kind, see again |7]. Regularity of h with 
respect to o is obtained first on [0, Q applying the Implicit Functions Theorem 
in the integral formulation, and then on [(, ( Q ] by classical regularity of Cauchy 
solutions with respect to initial conditions and parameters. □ 
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For the sake of clarity we will denote by v = v a this unique regular solution 
of (5.3) and will only consider a > 0. 



Proposition 5.2. For a > only three scenarios are possible for (5.3) at 
infinity: v(+oo) = +oo, v(+oo) = and v(+oo) = — oo. Moreover, if ( a = 

m — 2 

f m zs the first time where b(Q = cr, the following holds: 

1. If there exists (i > ( a such that f (Ci) > and f'(Ci) > ; then v(() > 
v(d) > on [d, +oo[ and v(+oo) = +oo. 

2. If there exists (i > ( a such that f (Ci) < and f'(Ci) 0, then v(() < 
v((i) < on [(i, +oo[ and t>(+oo) = — oo. 



Proof. When £ = +oo (5.3) reads —v" + B(^™^v = and therefore satisfies 
the classical Maximum Principle (B = est > 0). As a consequence either 
v(+oo) = ±oo, either v (+oo) = est, and the only possible finite limit is clearly 



f(+oo) = 0. Recasting (5.3) as [( m - 2 v']' = \ v we have by definition 

Q m — 2 

b(Q — a > for ( > ( a , and the rest of the statement is easily obtained 
integrating from (j to +oo. □ 

Proposition 5.3. For o > small enough v > holds on [0, +oo[ ; and 
v(+oo) = +oo. 



Proof. Denoting by v the solution for a = 0, (5.3) reads [("^v']' = B( m - 2 v. 
Using boundary conditions v(0) = 1 and u'(0) = it is easy to integrate from 
C = to C > and see that v(+oc) = +oo. This easily extends by continuity 



m-2 



to cr > small enough remarking that £ CT := (§) m — > + when a — > + and 



using the first case in Proposition 5.2 □ 



Proposition 5.4. For a large enough there exists a time = C-(°") > such 
that v(C-) < 0. 



Proof. Scaling (5.3) as t = cr( and y(t) = v f^J yields 

1 



ty + ay +y = /3t™-iy, y(0) = 1, y(0) 



a 



with y — %, a = ^2 G]0, 1[ and /3 = — ^=r - When /3 = we may look 

for the solution in the form of a power series, and an easy but technical 
computation shows that this power series takes negative values y a (t = 2a) < 
for any fixed a g]0, 1[. This easily extends by continuity to (5 > small, i-e to 
a large. □ 

We may now prove our main statement: 
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Proof, (of Theorem [2]) . By Propositions 5.3 and 5.4 the quantity 



a =: sup (a > 0, o' G [0, cx] =S> v a \) > o) (5.5) 

is finite and positive. We claim that cr is indeed the desired principal eigen- 
value: denoting by Vq the corresponding solution for a = a we have by con- 
tinuity that vq = lim v a > 0, and actually v > (unless t>o = 0, which 



'0 



contradicts t>o(0) = 1). Stability t>o(+oo) = is a consequence of Propo- 
sition 5.2: the scenario t>o(+oo) = +00 is impossible since otherwise there 
would exist a right neighbourhood of o"o in which v a > 0, thus contradicting the 

m — 2 

definition of cr . In order to obtain monotonicity, let Co ; — \b) m '■ by Propo- 
sition 5.2 v > must satisfy v' Q < on [£ ,+oo[ (otherwise t> (+oo) = +00, 
which we proved to be impossible). On [0, Co] monotonicity is a consequence of 
v 'o( Q ) = ~( m ~ 2)^0 < and Lv = (cx - b)v < 0, where L = -(-^ - 
is elliptic and has no zero-th order coefficient. □ 

5.2 Isolated eigenvalue 

We prove here that there are no other eigenvalues close to o-q. More precisely 
we show that a G [0,cto[=^ y(+oo) = +00 and that for a > o"o close enough 
v(+oo) = —00 holds. The key point is the following monotonicity: 

Proposition 5.5. For a G [0,o"o], the mapping a \-t v ff (.) is strictly pointwise 
decreasing on ]0, +00 [. 

Proof. Remark that ( = is irrelevant since we normalized v a (0) = 1 for all a . 
The proof is quite standard: if 0\ > a 2 G [0, a ] then t> 1; t> 2 > 0, a := ^ > is 
well defined and satisfies 

L[a] := -(v 2 a" - [ 2(v' 2 H ^— ) a = (oi - a 2 )v 1 > 0. 

\ m — 2 J 

Computing a(0) = 1 and o/(0) = (m — 2)(a 2 — a\) < shows that a < 1 at 
least in the neighbourhood of ( = 0. If ct(C) ^ 1 for some C* > then a would 
attain a non-negative local minimum point ( m g]0, £*[, which would contradict 
the classical Minimum Principle. Thus a < 1 <^ V\ < v 2 on ]0, +00 [. □ 

Proposition 5.6. For a G [0,ao[ we have that v a (+oo) = +00. 

Proof. For such a we have by monotonicity v a > v ao > hence by Propo- 
sition 5.2 either Vcr(+oo) = +00, either v a (+oo) = 0. Assume by contra- 
diction that ^(-l-oo) = for some a G]0,<7o[: both v a and v ao then satisfy 
the same asymptotic equation when ( — > +00 and thus behave similarly. 
Slowly increasing a from zero, we see that av a < v ao holds until a first 
critical value a = a* G]0, 1[ for which there exists a contact point ( m be- 
tween a*v a < v aa and v ff0 . On the other hand z := a*v a — v uo < satisfies 
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Lz = (cr - cr ) a*v a + d z < with L = -(-^ - + 6(C), thus con- 



<o 



<0 



tradicting the Minimum Principle. □ 

Remark 5.2. In the proof above we only used the fact that cr is associated 
with a (stable) positive eigenfunction v ao . As a classical byproduct we retrieve 



uniqueness of the principal eigenvalue cr : by Proposition 5.6 there are no 
eigenvalues at all for a < a , and if there existed an other principal eigenvalue 
o\ > o"o associated with a positive eigenfunction v ai we could repeat exactly the 
same argument and conclude that v a > v ai and w cr (+oo) = +00 for all a < o\, 
which fails for o = cr < °\ ■ 

Lemma 5.1. If a > o-q then v a > cannot hold on R + . 

Proof. Assume by contradiction that v a > for some a > a . Once again we 



must have either f cr (+oo) = 0, either v a (+oo) = +00. Owing to Remark 5.2 
above the first case is impossible, and therefore v a (+oo) = +00. Defining 
a := ^f- we have a > 0, a'(0) = (m — 2) (do — a) < and a(+oo) = +00. As a 

consequence a attains a (non-negative) minimum point at some ( m > 0, thus 
contradicting the classical Minimum Principle and 

L[a] := -(v a" - [2(v' + ) a ' = (a - a )v > 0. 

\ m — 2 J 

□ 

Proposition 5.7. For a > a close enough there holds v a (+oo) = —00. 



Proof. By Proposition |5.2 it is enough to show that for such a there exists a 



time C > Co- = (§) m such that v a (Q < and v' a (Q < 0. For any a > a we 



normalized v a (0) = 1 > 0, and by Lemma 5.1 v a (.) takes negative values at 
least somewhere in M + : the first time Co( cr ) > where v a vanishes is hence well 
defined. Since v ao > we must have by continuity Co(° r ) — > +00 when o — > o~q , 

m-2 

and in particular Co(°") > Co- since Co- — > (^) m < +00. By definition of Co 
we also have v a > on [0, Co[ an d ^o-(Co) = hence t^(Co) < 0, and actually 
v a(Co) < 0- Thus the desired behaviour for times C ~ Co(o") + and a > o"o 
close. □ 



5.3 Analyticity 



For cr > equation (5.3) has a left branch of solutions Vi(a, .) which are C 1 
at C = and normalized as Vi(a, 0) = 1 (see Proposition 5.1). According 



to Remark 5.1 there also exists for any a a right branch of stable solutions 
v r (o~, +00) = 0, which we normalize as v r (o~, Cq) = 1 for some Co > 0. In this 
section we show that these two branches are analytical in a, which we will 
need later for technical reasons. 
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Proposition 5.8. For any fixed Co > there exists R(Co) > such that 
a i — y vi(o~, .) is analytical into C 1 ([0,Co]) with radius of convergence at least 



Proof. If no singularity occurred at C = in (5.3) this would be a classical 
consequence of the linear dependence of the equation on a. We construct 
below a solution v G C 1 ([0, Co]) °f (5.3) which is analytical in a and such that 
v(0) = 1. By uniqueness in Proposition 5.1 this solution will agree with Vi(a, .). 
For a w o~q let us enlighten <r — o"o and rewrite (5.3) in the form 

Lv := -(v ' - m _ 2 +(b- cr )v = (a - a )v. (5.6) 

Looking for solutions in power series v (() = Xln>o( (T — cr o) n ' y n(C) leads to 

n = : Lvo = 
n > 1 : Lt> n = 

We naturally set t>o := v ao to be the principal eigenfunction in Theorem [2j 
and also require v n to be C 1 at ^ = and satisfy v n (0) = for n > 1 (so that 
u(0) = 1 in the end). For any g G C([0,Co]) it is easy to check that Duhamel 
formula 

/(C) = -v (0 [ - 1 [ I v^t^-'g^dt | da. 

yields the unique solution Lf = g satisfying the above boundary conditions, 
and that H/Hc 1 < C| Is'l |oo f° r some C(Co) > 0. Equation Lv n = v n ^% is 
therefore uniquely solvable as v n = L~ 1 v n _i and an induction argument im- 
mediately yields ||w n || c i < C n ||t>o||oo, ^ nus convergence in C 1 with radius at 
least R((o) := g^y. □ 

We have a similar result for the stable right branch: 

Proposition 5.9. For any fixed Co > there exists -R(Co) > such that 
o i — y v r (a, .) is analytical into C^([Co, +oo[) with radius of convergence at least 
fl(Co). 

Proof. The argument is almost identical, the only difference being that we want 
now (v n ) n >i to decay at infinity (whereas we had to handle the singularity at 
C = for the left branch). The corresponding Duhamel Formula for Lf = g 
reads now 



Co 




/(C) = MO I ~ : I / voit^-'g^dt | da, 

VQ{S)S m - 2 



and the rest of the proof is straightforward. □ 
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5.4 Functional setting multiplicities 



Let us remind at this point equation (5.3) 

u' 

C > : -(u" h bu = an, 

m — 2 

which is an eigenvalue problem for L := — — m^'ck ^ ^(0- Boundary 
conditions are (i) C 1 regularity at ( = and (ii) decay v(+oo) = 0. So far 
we did not define any clear functional setting, and this is the purpose of this 
section. 

Since v(+oo) = the largest possible space on which L could act is clearly 
a subspace of 

E:={feC([0,+oc]), /(+oo) = 0}, 

and our functional setting should also take into account the singularity £ = 0. 

Proposition 5.10. For any f G E there exists a unique solution v G E of 
Lv = f having C 1 regularity at ( = 0. This solution can be computed as 



+oo Is \ 

v(C) = L- 1 [f](C) :=*(C) / _ 2/ ^ -j- / tJWt^-VWdt d S , 



(5.7) 



where v(() > is the solution in Proposition 5.1 for a = 0. VFe /iawe moreover 
i/(0) = -(m-2)/(0). 

Proof. The C 1 regularity at = is easy to retrieve, and one readily checks 
that L -1 maps E to E (or in other words that t> decays at infinity if / does). □ 

This defines very naturally the domain 

D{L) :=L-\E) C E, 



where L~ 1 is given by Proposition 5.10 Any u G D(L) therefore satisfies both 
boundary conditions, and let us point out that if / = L~ 1 g with g G E then 
/ is of course a classical C 2 (]0, +oo[) solution of the corresponding ODE. 

Proposition 5.11. o"o is an eigenvalue of L : D(L) C E — )■ E with geometric 
multiplicity m 9 (o~o) = 1- 

This is of course consistent with the one-dimensionality in Theorem [1} 

Proof. For a = <7o the principal eigenfunction vq = v ao constructed in Theorem 
[2] solves the ODE, has the required regularity and decays at infinity, thus 
belongs to D(L) and is indeed a proper eigenvalue. Regarding it's geometric 
multiplicity we recall from Remark 5.1 that for any a there always exists 
a solution of the ODE which blows (exponentially) at infinity and therefore 
cannot belong to E. □ 
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The algebraic multiplicity will also be important in Section |6j 



Proposition 5.12. oo has finite algebraic multiplicity m a (o"o) = 1. 

Proof. If Vq denotes again our principal eigenfunction and Lq := L — o~q we 
need to show that Ker(Lp) = Ker(Lo) is generated by Vq; in other words, that 
there is no solution v G D(L) of Lv = Vo — (V — ^2 + (P ~ a o) v — v o- 
Assuming by contradiction that v is such a solution, Variation of Constants 
v(Q = «(C)' y o(C) an d boundary conditions lead to 



Decay vq(+oo) = (at least exponentially) implies that the integral above 
is absolutely convergent at infinity, hence that a'(C) ~ ° i —■ —oo 

with C = Jq +0 ° s m - 2 ~ 1 vl(s) ds > (let us also recall that vq > 0). Using the 
equation for vq it is finally easy to see that a' blows so fast at infinity that 
v = avo does too, thus contradicting v(+oo) = and v G D(L). □ 

Proposition 5.13. a = a G R an isolated eigenvalue in C. 

Proof. Rewriting Lv = av ^ —(v" — —^v' + b(Qv = ov in the divergence 
form 

/. 1 a' HO ° 



1 .1 1 

m-2 



we see that the problem is self-adjoint and it is enough to prove that ctq is 
isolated in IR. Let us remind that we are looking for eigenfunctions v G D(L), 
hence C 1 at ( = and decaying at infinity. As a result any eigenfunction v a 



must agree with the unique regular solution of Proposition 5.1 which according 



to Proposition 5.6 and Proposition |5.7 cannot decay if o is close to o~ . □ 



6 Asymptotic matching 

In this section we match the physical setting e > 0, k < +oo with the (formal) 
asymptotic problem e = 0, k = +oo of Section |5j 



For s of order k 1 m ~ 1 we constructed in Section [4] the solution of (2.7a) 



in £ = x/e coordinates with maximal decay at —00, normalized as f(£ £ ) = 1. 
This was done for e > and k < +00 (small and large enough) so we may 
indistinctly switch from £ to x coordinates. This left solution extends to x > x e 
and will be denoted below by ui(x) = ui(e, k, s, x). 

On the right side x = +00 we have p = 1, p' = p" = and the characteristic 



equation associated with (2.7a) reads — r 2 + c £ r + (k 2 — s — G'(l)) = 0, yielding 
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two exponential rates r ± . In the regime (|2.9|) with s of order k 1 <C k z 
(and c e — > Cq > 0) these are 



r ~ ±/c, 

e.k 



and there always exists a unique branch of stable solutions m(+oo) = which 
we will denote by u r (x) = u r (e, k, s, x). 

Instead of parameters k — > +00, s = O (k l ~™^^ — > +00 and x coordinates 

we will rather use 

1 5 

y Q, a = s-. — = 0(1), ( = x/5, 



k l ~^^ (m — 2)c, 



and (5.2) reads 

Lv := -qv" + (c - 2^) v' + ({q)^ - (m - 2)ca - q" - 5G'{6q)) v = 0, 

c = c £ >0, q = q £tS (C) = Ipe(SC). 

(6.1) 

In terms of these new parameters (e, 5) the frequency regime (2.9) becomes 
u := \e, S > : < £ " n ^ + ^ < 5 < 5 \, (6.2) 



where rj > is again exactly the one in Theorem [T], and we are now interested 
in the double limit (e, 5) — >■ (0, 0). 

Throughout this entire section we will write Vi(e, 5, a, () for the left branch 
coming from the cold zone and v r (e, 5, a, C) for the stable right branch normal- 
ized as tv (Co) — 1 (for some Co > to be chosen later). Those are simply u\ and 
u r expressed in ( variables and in terms of (e, 5) instead of (e, k). We will also 
denote by Vio(a, C) the regular solution of the formal asymptotic problem (5.3) 
(defined in Proposition 5.1), and v r o(a, () its stable right branch normalized 
as f r o(o", Co) = 1 (see Remark 5.1). We will prove that 

vi(e, 5, a, .) vio(cr, .) v r (e, 5, a, .) v r0 {<r, .) 

when (e, 5) — > (0, 0) in u. This should be no surprise since the asymptotic 
problem was precisely obtained in the (formal) limit (e, k) = (0, +00) or equiv- 
alently (e, 5) = (0,0). 

Our first step is to establish convergence of the scaled wave profile q e> g — > qo 
in the linear zone C £ [CejCsJj where Ce = x e /5 corresponds to the exit of the 
cold zone x = x £ = e 1_a and C<5 = xg/5 to the exit of the linear zone x = xg 
defined in (3.4). One can check that 

Ce — > and C<5 ~ +°° 
when (e, 5) A (0, 0), and the linear zone stretches to ]0, +oo[ in this scaling. 
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Proposition 6.1. We have convergence \\q Et s{C) 



[m 



2)c C|| 



the double limit (6.2) 



C 2 (Ke,( S ]) 



— > in 



sition 



The proof is technical but straightforward when taking advantage of Propo- 
3.1 and asymptotes (3.6). In order to retrieve the above C 2 convergence 



it is important that the exit of the cold zone was pushed far enough out of the 
boundary layer, thus allowing to neglect p"(x) in the linear zone. 



6.1 Regularity of the branches 

We start by establishing convergence of the left branches as claimed above: 

Proposition 6.2. Fix( > anda*: thenvi(e,5,a, .) — > vi (a*,.) mC 1 ([C £ ,Co]) 
when (e,5,a) — > (0,0, cr*) in the double limit (6.2). 

Note how we handle simultaneously the double limit (e, 5) — > (0, 0) and 
continuity with respect to a. 



Proof. For e, 5 > the left branch vi = vi(e, 5, a, () solves (6.1 ): 

Lvi = 0, L — -q^r 2 + (c - ~{m- 2)co - q" - 5G\5q) 

while the asymptotic left branch %i = ^(c^C) solves (5.3): 

d 2 d 

L v m = 0, L = -(m - 2)c C^ - c — + (m - 2)c (6(C) - o - *). 

Proposition |6.1| guarantees that the coefficients of L are very close to those of 
Lq when (e,5, a) — > (0,0,0"*). Boundary conditions (4.14) moreover read in ( 
coordinates 

Vl (Q = l, v' l (Q~-( m -2)a 

E,0 



with C £ — > 0. According to Proposition 5.1 the asymptotic left branch satisfies 

v lo (0) = l, v' l0 (0) = -{m - 2)a„ 

and we see that initial values are also very close. These two branches solve 
two very similar Cauchy problems, and should therefore stay very close on the 
interval [C e , Co]- The difficulty is here that the asymptotic Cauchy problem 
degenerates at ( = when e, 5 — > 0, thus preventing from applying classical 
regularity of Cauchy solutions with respect to parameters and initial values. 
In order to keep this paper in a reasonable length we only sketch the argument 
below. 

The proof relies on a stability argument for z := V\ — viq, which satisfies 
z(Q = o(l), z'(Ce) = o(l) and Lz = f with / = (L - L)v m = o(l). Using 
the very structure of the equation one shows that z stays small in C 1 norm for 
times ( G [Ce>C]> where ( > stays bounded away from zero when (e,5,cr) — > 
(0, 0, a*). This allows to step away from the singularity ( = 0, and we conclude 
on [(, Co] by classical regularity arguments. □ 
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We have a similar result for the cr-derivative: 



Proposition 6.3. Let z\ = and zio = For fixed Co > and cr* we have 
Zi(s,S,a,.) — > zjo(c*, .) in C 1 ([C e ,Co]) w/ien (e,8,a) — > (0,0,(7*) m the double 
limit (6.2). 

Proof. The proof is very similar to the previous one: differentiating Lvi = 
and LqViq = with respect to o we see that z\ and z\q solve again two very 
close equations, and Proposition 4.15 with t>2o(0) = 1 and v' l0 (0) = — (m — 2)cr 
show that initial values are also close in the above limit. □ 

The delicate issue for the right branches is to take into account the decay 
v r (+oo) = 0, which turns to be a slightly more difficult problem than the 
previous singularity at C = for the left branches. We first establish two 
decay estimates on v r , v' r at infinity and uniformly in e, 5, a: 

Lemma 6.1. There exist Co > and C > locally independent of a such that, 
if e,5 are small enough in the double limit (6.2), then < v r (e,5, o, Q < Ce^^ 
holds on [Co, +00 [. 

Proof. Let us recall that v r solves the elliptic equation Lv r = 0, with 



d2 , 



fc - A + [q^ - (m - 2)ca - q" - 5G'(5q) 



The key point is that, if Co is chosen large, the zero-th order coefficient 

Qo : = g^i _ (jn - 2)ca- q" - 5G'{5q) (6.3) 

O(l) 0(5) 

can be made larger on [Co,+oo[ than some positive constant locally inde- 
pendent of a: since q(Q ~ (m — 2)coC it i s indeed enough to take Co > 

m-2 

(3(m — 2)<r) m and thus ao(C) > (^ — 2)cocr > est > (locally in a). As 
a consequence we may apply the classical Minimum Principle Lv r = with 
v r{Co) — 1 and t> r (+oo) = to retrieve positivity. To obtain the upper esti- 
mate, remark that at infinity the asymptotic equation associated with Lv r = 
implies that v r behaves as e~ r< * for some positive r ^> 1: it is then easy to 

e,k,a 

check that v = Ce^^ is a supersolution for v r if C > is chosen large enough 
(but independent of e and S and again locally of a). □ 

We will also need a similar estimate for the first derivative: 

Lemma 6.2. There exist Co > and C > locally independent of a such that, 
if (£,5) are small enough in the double limit (6.2), then \v' r (Q\ < c 4 

gm-'2 

holds on [Co, +00 [. 
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Proof. Fix any Co > large enough as in Lemma 6.1 dividing (6.1 ) by q > 

we recast Lv r = as 



v" + av' r = bv r and 



q \ m—2 

b := i (q^> - (m - 2)ca - q" - 5G'(5q) 



and obtain the corresponding integral formulation 



+oo 



V'r(0 



exp [A(s) - A(()} b(s)v r (s)ds 



with A(s) := f s a. Taking advantage of the explicit form of a(C) combined 
with q(Q > q(Co) ~ (m — 2)cqCo an d q( s ) < = t it is easy to estimate 



exp [A(s) - A(Q] 



s > C > Co : < 



Uo 











J exp | 













< 



2 _ 5 



<1 



1 ^ O / 1 \ m-2 1 



{ \b{s)\<2q^- 1 <2{\) 



if Co > is large enough. Our statement finally follows from Lemma 6.1 
above. □ 

We have now for the right branch 

Proposition 6.4. For any a* there exits Co > large enough such that 
v r (e, 5, a, .) —> tvo(o"*, •) in C 1 ([Co, +oo[) when (e, 5, a) — > (0, 0, cr*) in the double 



limit (6.2); Co can be chosen locally independent of a*. 



Proof. From Lv r = and Lqv t q = we see that w := v r — t> r o satisfies as 
before Lw = f '■= (Lq — L)v t q on [Co, +oo[. Contrarily to the left branches the 
coefficients of L are not close to those of L on the whole interval [Co,+oo[, 



but only on [Co, Cs] (this is because of our C-scaling, see Proposition |6.1 ). From 
(5.3) it is however easy to see that v r0 , v' r0 , v" decay at least exponentially at 



infinity, and locally uniformly in a*. According to Lemma Chi and Lemma [672 



+oo it is 

e,5 



we know that v r ,v' r also decay exponentially, and since C<5 = ^ 

clearly enough to prove our statement on the sub interval / = [Co, Cs] • 

Remark that (i) the zero-th order coefficient of L can be made large and 



positive, see (6.3), (ii) C 2 convergence of the wave profile q £} § — > go implies that 
f = (Lo — L)v r o is small in L°° norm on / and (iii) w is exponentially small 
on the right boundary and vanishes on the left one, see Lemma 6^ Using 
w := |w|(C<5) + | l/l |oo = est as a supersolution Lw > f we obtain \w\ < w on 
/ = [Co? 0] an d w = o(l) by construction. As a consequence v r and f r0 are 
close in L°° on [Co, +oo[, and we are left to control the derivatives. 
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According to Lemma 



6.2 



we have that lull < 



c 



-Ct 



s m - 



o(l) on [C5,+oo[ 



and v' r0 also decays exponentially: thus \w'\ = o (1) uniformly on [Q, +oo[. In 



e,s 



order to retrieve the same estimate on [Co, Cs] we integrate the equation for w 
backward in time from £5, where w, w' are small. □ 

A similar result holds for the a derivative: 



^ and z r o = for any fixed cr* t/iere exists 

) — >■ z r o(c*, •) wi C 1 ([Co>+oo[) w/ien (e, 5, a) — ► 



Proposition 6.5. Lei z r = 

Co > sfic/i £/ia£ z r (e,5, a 
(0,0,a*) m i/ie double limit (6.2); Co can be chosen locally independent of a 



Proof. The proof of Proposition |6.4| is easily adapted, the first step being to 
prove two decay estimates |-z r (C)l ^ Ce~^ and |^(C)| < ^ e ~^ 011 z = ^ f° r 
some fixed exponent a = a(m) > (similar to Lemmas 6.1 and 6.2 above). □ 



6.2 Evans function and construction of the eigenvalue 

In this section we construct the principal eigenfunction a = o~(e, 5) connecting 
maximal decay 00) = and v r (+oo) = in the frequency regime (6.2). 
Rewriting (6.1) as a first order system 

dY „ „ f v 



we denote by Yi(e, 5, a, .) and Y r (e, 5, a, .) the previous left and right branches 
and define the Evans function 

E(e, 5, a) := det [Y^e, 6, a, Co), Y r (e, 6, a, Co)] (6.4) 

for some Co > to be chosen later. Clearly the connection holds if and only 
if E(s,5,o~) = 0, thus reducing the problem to finding the zeros (e,5,a) of 
E. Let us recall that the frequency open set u is defined in (6.2), and let us 
introduce 

to := u U {(0,0)}, Q := cox]a - A,a + A[, Cl := iux]cr - A, a + A[ 

(6.5) 

for some A > 0. For (e, 5, a) G Q (hence e, 5 > 0) there are no singularities in 
the equations: the two branches Yj, Y r are well defined and regular in (e, 5, a), 
and the Evans function E(e, 5, a) is therefore well defined and C 1 on Q. The 
set Q however contains critical points of the form (0,0,0"*), corresponding of 
course to the asymptotic problem e = 0, k = +00. 

Proposition 6.6. Fix Co > large enough: the Evans function extends to 
a continuous function on Cl, again denoted by E(e,5,a). This extension is 
continuously differentiable with respect to a in Cl. 
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Proof. Since we consider bounded o G]cr — A, ctq + A[ Propositions 6.2 and 



6.4 hold for some Co > large enough but independent of e, 5, a (let us stress 



that all the results in Section 6.1 were locally independent of a). The branches 

Yi,Y r thus continuously extend to Yi Q ,Y r0 when (e,5,a) — > (0,0, er*), and this 
obviously extends E to Cl setting 

E(0, 0, a) := det (Y m (a, Co), Y r0 {a, Co)) • 



By Propositions 5.8 



and 



and Propositions 6.3 and 



5.9 



6.5 



the mapping a i— > .E(0,0, cr) is analytical in cr, 
finally imply continuity of the cr-derivative with 



respect to all three arguments (e, 8, a) at any asymptotic point (0,0,0"*). □ 



According to Theorem [2] and by definition of the extension E(0, 0, .) we 
have of course 

£(0,0, oo) = 0; 

in order to apply an Implicit Functions Theorem we need to check the usual 
condition: 

Proposition 6.7. E satisfies §f (0, 0, cr ) ^ 0. 

Proof. The argument is directly inspired from [T], Lemma 6.2 p. 194. 
5 = the first order system corresponding to (5.3) is 



For e = 

dY 
1( 



A(a,0Y, Y 



v 
v 1 



A(a, C) 





b-o 

c 



1 



1 



(m-2)C 



and the Wronskian W(o~, () = det (Yi (a, C), Y r0 (a, C)) satisfies 

C^W(a,() = est = Cr 2 W(a,Co) = C ^^(0, 0, a). (6.6) 

6(C) as an un- 



In the light of Section 



5.4 



we consider L 



1 d 

' dC 2 rn-2 d( 

bounded operator L : D{L) C E — > E, of which <7o is a (principal) eigenvalue. 
For cr 7^ o"o and / £ E computing the resolvent R(cr; L)f amounts to solving 
— a)v = f with boundary conditions v G D(L) corresponding 
to C 1 regularity at C = and decay v(+oo) = 0. Using Variation of Constants 

mputed as 

avn + fivro), (6.7) 



and (6.6) this unique solution can be computed as 

1 



with 

«[/](C) 
PifKO 



t^=2<f>(t)v' r0 (t)dt, 
t^${t)v' l0 (t)dt, 



and $(C) := 



t^2~ 1 f(t)dt. 
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According to (6.6) and (6.7) the resolvent therefore reads 



V/e£, R(a;L)f = v 



E(0,0,a)Co 



—M{a)f, 



(6.8) 



where M(o~)f 
tions 



(aviQ + /3v r o) and a [/],/?[/] are defined above. By Proposi- 
5.8 and 5.9 the mappings a i-» Vio(a, .) and o \-t v r0 (cr, .) are analytical, 



and the family of operators M{a) is therefore analytical. We constructed o"o 
so that the principal eigenfunction t> = fzo(co; •) = v r o(ao, .) G -D(L) (up to 
some proportionality factors): it is then the easy to compute M(<jq)vq ^ 0, 
and M(o"o) is consequently non-trivial. 



By Proposition 5.13 the principal eigenvalue is isolated in C, and the resol- 
vent R(cr; L) is thus classically meromorphic in some complex neighbourhood 



of o"o. Since E(0,0,a) and M(cr) are analytical and M(ctq) ^ 0, (6.8) shows 
that the order of the pole a = (Tq in R(a; L) equals the order of the zero 
£7(0, 0, o"o) : by classical functional calculus arguments the order of the pole 
in the resolvent is exactly the algebraic multiplicity (see e.g. |12j). which is 
m a(°"o) — 1 by Proposition 5.12 □ 



We may consequently apply the following Implicit Functions Theorem: 



Theorem 3. Lettu be as in ( 6.2[ ): if 5$ > is chosen small enough there exist 
B g]0, A[ and a function a(e, 5) defined on u such that 



(e, 5, a) G ux]o~q — B, a + B\ 
E(e,5,(t) = 



(e,5) G u 
a = a(e, 5) 



2. a is C 1 on u and a(e, 5) — > ctq when (e, 5) — > (0, 0) in u. 



Note that, due to the algebraic form (6.2) of u>, taking 5 < So small enough 
is just a vicinity condition (e,5) ~ (0,0) in u, that is to say u (e,k) G U and 
k > k large enough" as stated in Theorem [TJ 

Proof. This is a non-standard version of the Implicit Functions Theorem since 
E(e, 5, a) is not C 1 in some open set, which Cl is not. The argument is based 



on the usual iterative fixed point and strongly relies on Proposition 6.7, but 
will not be detailed here. □ 



6.3 Proof of Theorem Q] 



We just defined the function cr(e,5) on u so that E(e,5,a(e,S)) 
other words so that the connection v(±oo 



0, or in 
a(e,8). 

We denote below by u e ,a(C) the corresponding eigenfunction and by u £ ^{x) its 



holds in (6.1 ) if a 



alter-ego in x coordinates, which is a solution of (2.7a) for e > and k < oo. 
Our previous scaling a = s 1 t and 5 = \ \ corresponds now to 

m — 1 



(m— 2)c e k 



s(e, A;) := (m — 2)c £ cr ( e 



_ i 

m— 1 
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which is of course the principal eigenvalue in Theorem [T] This immediately 
yields the desired asymptotic relation s ~ 7o/c 1_m - 1 with 70 := (m — 2)co<To. 

e,k 

For a 



a(e,5) ~ o"o > the asymptotic equation of (6.1) at +00 leads to 
an unstable characteristic exponent r + > and the stable subspace v(+oo) = 
has therefore dimension 1. This shows that the eigenspace for s = s(e,k) is 
1-dimensional as stated in Theorem [T] We may have equally argued that the 
space of maximal decay solutions at —00 is also 1-dimensional for s = s(e, k) ~ 
7oA: 1_m - 1 , see Section [4] and in particular the discussion about B w spaces page 

EU 

We prove below the claimed positivity u £ ^{x) > <^ Ve,s{() > separately 
on three intervals: in ( coordinates on [£ £ , Co] (for some Co to be chosen later), 
again in ( coordinates on [Co, +00 [, and finally in x coordinates on ] — 00, x £ ]. 



C G [C Co] Since a(e, 8) ~ ctq Proposition 6.2 guarantees that v E> s(C) * s 



close to the asymptotic principal eigenfunction Viq((Tq, () uniformly on [( £ , (q\ 
for fixed Co > when (e, 5) — > (0, 0) in the double limit. By construction we 
had v m (a , () > on E + , thus v £jS (() > on [( £ , ( ]. 

For a = a(e,5) equation (6.1) is of the form L [v £j s] — 0, 



C e [Co, +00 



where L is uniformly elliptic. As already discussed the zero-th order coefficient 



a o(C) given by (6.3) can moreover be made positive on [Co,+oo[ if Co > is 
large enough: we conclude applying the classical Minimum Principle with 

Ve,s (Co) > and v £iS {+oo) = 0. 



x G — 00, x. 



Here the classical Minimum Principle is unfortunately not 
available at once since the zero-th order coefficient in Lv = is not positive, 
see again (6.3) with q{— 00) = e/5 1. In order to keep our notations 



light we omit below the subscripts and one should understand in the following 

v(() = v £jS ((), u(x) = w e ,fc(x), p = p £ {x), c = c £ and s = s(e,k). The 
asymptotic expansion in the cold zone v 
on ] — 00, x £ ] 



Vq + evi +EV2 reads in x coordinates 



u + EUi + EU 2 , 



u =p'/p'(x e ) > 0, 



and the leading order uq is positive. According to (4.9) the next orders eui, eu 2 
are small on 1 — 00, xJ, 



SE 



1-a 



< C lS £ 



1-a 



(6.9) 



\eu 2 \\l°° < C\\eu 2 \\ B 2 n < C 2 (s£ la ) 2 ' 

1 11 11 11 - LJ W ^Q \ ' 

and Uq > should suffice to ensure positivity of the whole expansion at least 



for times x < x £ . The asymptotic study of (2.7a) at —00 however shows that 



C e + 



U 



~ e 
-00 



< r 



c 2 £ + Ae ( k 2 E™- 2 



2e 



< 
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whereas the planar wave decays as p'(x) oc eE 1 C e r+x , and the lower order 

— oo 

terms eu\ + eu-i unfortunately dominate w when x — > — oo. We prove below 
that eu\ +SU2 stays negligible on some interval [xi, x e ], and use then a suitable 
comparison principle on ] — oo,Xi\. 



1. Fix A > G\ with G\ > as in (6.9): monotonicity u' = p"/p'(x £ ) > 



and se 1 a <C 1 define a unique time X\ G] — oo, x e [ such that 

u (xi) = Ase l ~ a , 



(6.10) 



and estimate (6.9) immediately implies 



x G [xi,x E ] 



u> ^uo^ -HewiHoo— Halloo 



> 0. 



2. Let ro := I s - > and := e rolE : using definition (6.10) of xi and 



asymptotes (3.3)-(3.5) we see that p ~ e, p' = o(l) and p" < Cse a hold 
on ] — oo, Xx] and hence 



-pr + c 



2p' 



+ ( k 2 p 



P 



3 r x 



> C [-er 2 + c r - s£- a ] 



a | ^rox 



> c 



4s 



— se 



> (s^" < 1) 



The function w := u/<& satisfies an elliptic equation Lw =: Lu = 0, 
where L has thus positive zero-th order coefficient ao=£[ < ! ) ] > 0. Since 
u ~ e r+x decays faster that $ = e r ° x (r + ~ cq/s > r ) we have 

— oo 

w(— oo) = 0, and also w(xx) = u(xi)/$(xi) > 0: we conclude apply- 
ing the classical Minimum Principle on ] — oo,X\]. 

In order to fully establish Theorem [l] we only have left to prove that 

(i) s(e, k) is the smallest eigenvalue (not necessarily principal), 

(ii) there exists no other principal eigenvalue. 

The proof of (i) is very classical and relies on the positivity of the prin- 
cipal eigenfunction u and the fact that, the smaller s, the faster a potential 
eigenfunction must decay on both sides (this can be seen looking at the char- 
acteristic equations at x — ±oo). Assume by contradiction that there exists 
a non-trivial eigenfunction u associated with some s < s(e, k): then w := u/u 
is well-defined, satisfies an elliptic equation Lw = with zero-th order coef- 
ficient (s — s)u > 0, and decays on both sides w(±oo) = 0. The classical 
Maximum/Minimum Principles finally imply w = 0. 

The proof of (ii) relies on (i): if there existed an other principal eigenvalue 
s > s(e, k) associated with a positive eigenfunction u we could apply (i) with 
s instead of s and conclude that s is the smallest eigenvalue. 



31 



References 



[1] J. Alexander, R. Gardner, and C. Jones, A topological invariant 
arising in the stability analysis of travelling waves, J. Reine Angew. Math., 
410 (1990), pp. 167-212. 

[2] H. Berestycki and B. Larrouturou, Quelques aspects mathema- 
tiques de la propagation des flammes premelangees, in Nonlinear partial 
differential equations and their applications. College de France Seminar, 
Vol. X (Paris, 1987-1988), vol. 220 of Pitman Res. Notes Math. Ser., 
Longman Sci. Tech., 1991, pp. 65-129. 

[3] H. BERESTYCKI AND L. NlRENBERG, On the method of moving planes 
and the sliding method, Bulletin of the Brazilian Mathematical Society, 
22 (1991), pp. 1-37. 

[4] E. Beretta, J. Hulshof, AND L. A. Peletier, On an ODE from 
forced coating flow, J. Differential Equations, 130 (1996), pp. 247-265. 

[5] L. A. Caffarelli, J. L. Vazquez, and N. I. Wolanski, Lips- 
chitz continuity of solutions and interfaces of the N -dimensional porous 
medium equation, Indiana Univ. Math. J., 36 (1987), pp. 373-401. 

[6] P. CLAVIN, L. MASSE, AND J.-M. ROQUEJOFFRE, Relaxation to equi- 
librium in diffusive-thermal models with a strongly varying length-scale, 
Communications in Mathematical Sciences, 9 (2011), pp. 127-141. 

[7] E. A. Coddington AND N. Levinson, Theory of ordinary differential 
equations, McGraw-Hill Book Company, Inc., New York- Toronto-London, 
1955. 

[8] L. MASSE AND P. CLAVIN, Instabilities of ablation fronts in inertial 
confinement fusion: A comparison with flames, Physics of Plasmas, 11 
(2004). 

[9] A. MELLET, J. NOLEN, J.-M. ROQUEJOFFRE, AND L. Ryzhik, Stability 
of generalized transition fronts, Comm. Partial Differential Equations, 34 
(2009), pp. 521-552. 

[10] J.-M. ROQUEJOFFRE, Eventual monotonicity and convergence to travel- 
ling fronts for the solutions of parabolic equations in cylinders, Ann. Inst. 
H. Poincare Anal. Non Lineaire, 14 (1997), pp. 499-552. 

[11] J. SANZ, L. MASSE, AND P. CLAVIN, The linear darrieus-landau 
and rayleigh-taylor instabilities in inertial confinement fusion revisited, 
Physics of plasmas, 13 (2006), p. 102702. 

[12] A. E. Taylor and D. C. Lay, Introduction to functional analysis, John 
Wiley & Sons, New York-Chichester-Brisbane, second ed., 1980. 



32 



